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GENERALIZED INDICIAL FORCES ON DEFORMING 
RECTANGULAR WINGS IN SUPERSONIC FLIGHT! 


By Hanvarp Lomax. Frinaten Bo FUL Le and Lowa Sit onK 


SUMMARY 


aLomethed is presented fur determining the tmne-de pendent 
Jlow overa rectangular wing moving with a supersonie forward 
speed and undergoing small rertieal distortions expressible as 
polynomials rneolring spanwise and chordwise distances. The 
solution for the velocity potential ts presented in a form analogous 
to that for steady supersonic flow haring the familiar reflected 
aren” concept diseorered by Evrard. Particular attention ts 
paed to mdicial-lype motions and results are expressed in terms 
Numerteal results for Mach 
numbers equal to Pt aud t.2 are given Jor polynomials of the 
Hirst and filth degree tn the chordwise and spanwise directions, 
respeeticely, ona wing having an aspeel ratio of 4. 


of generalized indicial forces. 


INTRODUCTION 


One of the basie problems arising in the analysis of wing 
flutter boundaries is the calculation of the aerodynamic forces 
on wings undergoing small but arbitrary spanwise and chord- 
wise distortions. When the wing aspect ritio is large (actually, 
when the distance between spanwise nodal lines is large), 
these forces are usually estimated by some strip theory in 
whieh the loading on each spanwise section is approximated 
from that on a two-dimensional wing having the same chord- 
wise distortion. This report is concerned with low-aspect- 
ratio reetangular wings for which tip effects are important 
and the full three-dimensional theory must be used. 

The exact linearized solution for the forees on thin: reetan- 
gular wings (limited, however, to the range where effective 
aspect ratio (Q vf 1 A) is 2b) traveling at) supersonic 
specds has been presented by both Gardner (ref. 1) and 
Miles (refs. 2 and 3) in terms of multiple integrals involving 
arbitrary surface undulations. However, the use of such 
solutions in evaluating. numerically say, the forees induced 
by specitie wing distortions still presents some difficulties. 
It is the purpose of this report to discuss certain techniques 
that ean simplify the labor involved in these caleulations and 
to present numerical tables for the forces induced by a class 
of surface deformations, a class general enough to represent 
the first few mode shapes of rectangular plates. 

Mathematically the problem is to find and analyze a solu- 
tion to the four-dimensional wave equation 

1 
Cro Fy E22 2¢4°=9 (la) 

do 
(where dy is the speed of sound. ( is the time, and 7.4.2 are 
space coordinates) that satisfies the appropriate boundary 


The particular form of the solution toa be 
analyzed differs from those presented by Gardner and Miles 
but its development ts based on the method due to Gardner. 

Hadamard (ref, 4: studied a generalized form of equation 
(lad tn which the number of dimeusions was arbitrary. His 
solutions to these generalized equations are fundamentals 
different, depending on whether the total number of dimen- 
sions is odd or even. En fact. the methods Hadamard de- 
veloped apply directly only to equations for which the total 
number of dimensions is odd. Solutions for the even cases 
(such as eq. (lad) are determined by a “method of descent’: 
that is. the solution for the next higher odd-dimeusioned 
equation is fonnd and then reduced Ts asad ibeperle ad 
of) one dimension. ft is apparent, however, that sueh oa 
technique is in itself by no means unique. ‘Phus, Hadamard 
found the solution to equation (lay by deseending froma 
solution to the equation 


conditions. 


] 
ga en) (hh) 


Fir ~ Fy” Fs: a¥ 
i) 


but there are many other partial differential equations and 
groups of partial differential equations governing a five- 
dimensional (r.y.2.€.0) space all of which satisfy equation 
(in) ina phine &- constant. Gardner discovered a set of 
equations containing equation (lay ina & - constant plane 
which are simpler than equation (lai in that solutions could 
be found and adapted to the boundary conditions for time- 
dependent motion by methods well known to aerodynamicists 
who have studied the flow about wings in steady supersonic 
flight. This is the essential part of Gardner's contribution 
and it represents the technique upon whieh the development 
of the solution presented in this report is based. Aetually, 
Gardner first applied a Lorentz transformation to equation 
(lay and then used his method outlined above. ‘The appli- 
eation of sucha transformation is unnecessary and fins the 
disadvantage that the resulting coordinates have lost. their 
direct’ physical signifieanes. We will apply Gardner's 
method of descent directly to equation (lay and then proceed 
to analyze the solutions so obtained, 

In order to simplify the analysis as much as possible, we 
will limit solutions to the plane of the wing, and, further, 
consider only indieial-type boundary conditions: in other 
words, unsteady motions in whieh the wing attains instan- 
taneously, at the time zero, a certain spanwise and chord- 
wise distortion which is thereafter fived. [Tt is well known 
that the transient responses to these indicial motions ean be 
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Used, ina superposition integral, to obtain responses to mus 
other types of unsteads motion: it partienlar, responses to 
the harmonic oscillations of nonrigid wings. 

Finally, the principal interpretation of the results will be 
made in terms of venevalized forces, since these can be used 
directly in either Mutter or gust studies. and it will be shown 
that the amount of labor required to calculate such forces is 
reduced by using reciprocity relations derived frou: the | 
general theorems presented in reference 5. 


LIST OF IMPORTANT SYMBOLS 
a aspect ratio 
ha speed of sound 
in amplitude of indicial-downwash distribution 


mee eq. (ais 
beta function cSee eq. (BlAa).) 
Boooetpo gq) incomplete beta function (See eq. (BISb).) 


Curia) influence function for effect of side edge (See 
eq EATON 
“ tift covtticient, fl 
or 
1. indicia! lift coefficient due to angle-of-attack 
change, without: pitching, ©), OU 
e Vaan { 
e./ indicial lift coefficient due to pitching for a wing | 
” 1 
rotating about its leading edge, C, hs | 
4 OG gen | 
O., pitching-moment coefficient, positive when trail- 
ing edge tends to sink relative to leading edge, | 
moment H 
qoNe i 
Cae indieial pitching-moment coefficient due to 
angle-of-attack change (without pitehing) ! 
Oey 
measured about the leading edge. Cn, = > : i 
Oo asu 
Ca! indicial pitehing-moment coefficient due to | 
pitching measured about the leading edge for 
au wing rotating about tts leading edge. Cg 
Ot 
~ OY yon 
e wing chord 
Preity generalized indicial force coefficient (See eq. | 
(3th) ; 
ran generalized: indicial force coefficient: (See eq. | 
(37 ' 
huryt distance of wing camber line from >=0 plane | 
M Mach number 
Ap | ‘ se a . 
oading coefficient (pressure on the lower surface 
“ Minus pressure on the upper surface divided | 
by free-stream advnamie pressure ) 
( i ) binominal coctlicient, ( i )= be 
um m m! (n— my! 
: : ye 8 l 
q dimensionless rate of pitching, UC 
; hee 
Yo free-stream dynamic pressure, 5 ole 
TP generalized coordinate 
Q. generalized foree corresponding to the gener- 


eralized coordinate q- 
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real part of 
ar sen Fa 
Vit tye 
vie ne By ome 
Wile semispan 
Wing area 
area of acoustic plan form 
aren of reflected acoustic plan form 
at 
time 
t 
4 
ae Mt 
3 
Wing Kinetic energy 
Wing potential energy 
forward speed of wing 


Ce. 


vertical velocity 

Cartesian coordinates, fixed relative to the fluid 
at infinity 

coordinates with origin on center of wing leading 
edge (See fig. 13.4 

coordinates with origin on center of wing leading 
edge at time zero (See fig. 14.) 

r 

‘ 

Mar-t 
B 

(4, —y bln n’) 

angle of attueck (angle between flight: path and 
plane of wing), radians 

vif ot 

wing angle of piteh relative to horizontal, posi- 
tive when trailing edge ties below leading 
edge, radians 

coordinate measuring fifth dimension 

free-stream densitv 

velocity potential 

portion of velocity potential induced by sources 
in acoustic plan form 

portion of velocity potential induced by presence 
of side edge 

potential funetion in five-dimensional space 

Subseripts 

regions inanr, & plane (See fig. 7.) 

upper side of wing, 2-04 

singularity (e. g., source) position 


V71T regions on wing shown in figure 4 


STATEMENT OF THE PROBLEM 
THE GOVERNING EQUATION 


Assuming a wing’s vertical motion is of such a nature 
that the velocities induced in the fluid are small relative to 
the magnitude of the wing's steady forward motion, the 
normalized form of equation (la) 


rr Say * Err See iu ile) 
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where @ ail can be used as the governing partial differential 
equation of the tlow field. 
tletermination of the velocity potential when the body or wing 
in question moves through the fluid, the axes remaining 
fixed with respeet to the sill thatd infinitely distant from the 
origin, For convenience we place the wing leading edge 
on the y axis at ¢) Oo and the side edge on the x axis, The 
wing thes atacconstant forward cin the negative v directions 


This equation applies to the 


speed so at subsequent times the feading edge ties aloug the 
fine wc Mit, where VZois the Mach uuamber, and the side 
edge move atoug thee anis as shown in figure f. 


Rietiee 2. Wine ti fined: .Gordinate: <3 stern, 


THE BOUNDARY CONDITIONS 


The (uid velocity normal to the surface of a solid: moving 
Ina frieGionless flutd must be zero. Hf the equation of the 
solid’s surface ts represented by 


Gory tod 
this boundary condition can be expressed mathematically, 
im terms of the coordinate system used in equation (le), as 


dU, dg dt 
or Ooror 


Oe OF ALA 


herr 
oy ay): Oo: 


+ 


Consider a thin surface near the 2-0 plane. Phe equation of 
the eamber line of this surface can then be expressed in the 
form 


Gayot ra htary ly 0 


and, assuming that thickness and lifting effects ean be 
separated linearly, the boundary condition for the camber 
line becomes 


Og Oh eee Oh O¢ 


oF 4 
ote. ort oydy Oz 


Ft te ean eal 


i 


Hf the derivatives of A with respect to eae of Che coordiiattes 
are small, the two middle terms can be neglected and the 
expression for the boundary condition reduces to 


Bf a 
We wish to simulate a rectaneular wing deformed ii- 
dicially by bending in the spanwise and chord wise directions, 
Fer this purpose. on the portion of the > 0 plane ocenpied 
by the wing phan form, the vertieal velocity. which deterniities 
the wing shape aecording to the previous equation. is assumed 
to have the form 


() io 
if; = = a. r ak ) ( : ) var 


Where eis chord length. a, is a constant and / and a are 
Integers > 0. 

The expression oe) Wt isused so that for/ Othe tangent 
to the wing camber tine at the leading edye is tangent to the 
eht-path angle of the leading edge. Consider, for example, 
the ease / lon 8 Phe dowuwash 


a 


w “ure Mh 


represents an infinite class of surface shapes tava 7 che Corn 


ia 


Avay ty tre Mtr farp| 2) 


2el “| - 
where fiurge isan arbitrary finetion and Ais. bx definition, 
the distaner of the wing’s camber Tine from the 2 plane. 
Since, within the accuracy of linearized theory. the solution 
for the flow about the wing depends only upon the value of 
aw tay), the londing on adl the wings represented by the 
thove equation is the same. 
Het us inspect the two special cases 


th) fi a yp tind 
i fing) 0 
For ease oi) 
dN E 


‘Del ” sine Ae 


Airy 


and the wing isa flat plate pitching at a uniform rate about 
its leading edge whieh is following the flight path 


_ tw 


Wine Hd 
alia 


as shown ? in figure 2. Henee. at time # the tangent to the 
fleht path of the leading edge is 


dW, dU al! 
sed i c 
The slope of the leading edge of the plate at the same time is 
he yl! 
or e 
and the two slopes are seen to be equivalent. 


The zs 
clear, A 


fe in both fleures 2 and 3 is purposely distorted inorder te make Che drawings 
sumption used insetting up the houndars value problern, by mes fowhich 
termined. wis that the surface of the wing must retain near the 2 =O plane, 
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Poorer 200 Flat plate pitching at noifern: rate about leading edge 


For ease ae 


Oe. Mite 


hisy.ts - : 
> el 
and the wing isa plate which obtained a sudden parnbatie 
camber at ¢. O, a shape tt maintained thereafter as shown * 
In fieure 3. 


x 
~ 
' 
we 
‘ 
“ 


Fieve 38, -Plate with parabotic camber. 


The problem is near, so it will be sufficient to determine 
a solution for arbitrary Zand a. and then add results for any 
combination of terms desired. Thus, the complete 
houndary conditions to be studied are 


ats 


y 
c 


wey dys (la) 


oF am r 


yey 
Jand, since the loading is zero over 
portion of the plane 


over the wing plan form 
the remaining 


off the wing (2b) 


sinee the loading is given by 


Ap 


4 70 
qu “I “tl ; 


SOLUTION FOR THE POTENTIAL 


Figure 4 shows the wing plan form on the surface of which | 
the potential is required, together with the system of axes: | 
also, traces in the 20 plane of the wave system set up by 
the indieial motion of the wing are indicated. The wave 
pattern for onty two edges is shown: the flight speed is super- 
soute so the trailing edge has no effect on the velocities in- 
duced over the wing surface, and the results are valid) fin 
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their entirety) only for gf > Eso the opposite edge either dias 
no effeet or one that ean be incorporated by simple super- 
position, 

The wave traces divide the wing area into several regions. 
indicated by the Roman numerals. in each of whieh dhe 
wmilytical formulation for the potentuel is different. Region 
7 consists of that part of the wing where the effeet of neither 
the side edge nor leading edge has vet been felt. 
77. the stde-edge influence ts acting cthe line yo Gis the trace 
of the starting evlindrieal wave from the side edywe v0 
Region /77 is the part within the 
sturting evlindrical wave from the leading edge, but outside 
the tatluence of the side edge. 


In region 


but not the lending edge. 


This region. and region ©. 
are further subdivided for reasons that will appear tater. 
Region 7s a compound region; potential there ean be 
found by adding the potentials for regions // and /7/ and 
subtracting the potential for region /. 
of the portion of the wing within the spherieal wave orig 
The flow over the pact of the 
wings comprising regions V/oand V/7 has reached a steady 
state relative toca pout ou the wing, and the potential there 
is just that for the corresponding parts of a rectangular wing 


Region Vo consists 


nating at the wing corner. 


with the proper downwash distribution in steady motion. 
Finally, region V/7// is againa composite region, its potential 
being the sunt of potentials for regions 7/7 and V// fess the 
potential for region V/, 

All the regions just: listed. with) the exception of region 
Voare actually governed by the three- (total) dimensional 
wave equation and the potential therein could be obtained 
by methods applicable to this simpler equation, However, 
in this report we shall present a unified approach and the 
problem will be solved by the same method in all regions. 


aoe 


fay e3ge so tote 


Fictke +. - Regions used in the analysis of a rectangular wink it 
supersoni¢ unsteady motion. 


REVIEW OF KIRCHHOFF'S FORMULA 


The solutions developed in’ the subsequent: sections are 
more clearly interpretable tf they are compared with certain 
known results that have already been determined for the 
indicial motion of nontifting wings with svomet rent! thickness 
distributions or lifting surfaces with all supersonie edges. 
The purpose of this section is simply to review briefly some 
of these latter results, 


GENERALIZED 


As in steady-state: wing theory, there is a formula for 
time-dependent flows that relates Che velocity potential to a 
distribution of time-dependent sources and doublets over a 
This formada is due to 
KRirchholf, and some of its aerodynamic uses are discussed 
Kirethot?’s result is immediately appheable 
inthe study of unsteady lifting-surface problems when the 


certain region in the wing plane. 
in reference 6. 
potential can be represented by sources alone. that is, when 


the tipper and tower surfaces of the wing do not taternet, as 
isthe ease tu revions J. 227, and Vioof tiga «4. 


Kirehhotf's formula for souree distributions can be 
Written 
bye des : 
Un, Ue {| | Edina (4) 
2a toor. 
mE 
where 
r ‘roo! (yy 


The brackets on a, indicate that the retarded value is to be 
taken 


wed oP) 


fur| 


und No indicates that the region of integration ts the acoustic 
plat form corresponding to the event crazier. 
veptsatre discussed at length in reference 6. 


These cone 


As has been pointed oat, equation (3+ holds for each of the 
revions 7. Eff and V7. but the area of integration S, differs 
considerubhy from one of these regions to another, Consider, 
for example, the determination of ¢ for region L77, denoted 
Part oof the boundary of the acoustic phar form S, 
is fond hy eliminating To between the equation of the 
leading eda. ay MET. and the expres-ion 


Sry. 


weary ype ait ut TA 
Which gives the outer boundary, at “time fof all the dis- 
turbanees that, operating at time’ 7) can produce an effect 


at the point cyt This boundary is the ellipse 


( - Fy ody )tuy- met, (da) 
where 
Mir -t ro Mt 
Pgh? hoe 
3 B 


If the point cay tes within the evlindrical wave from the 
leading edge, that is, feuect. the ellipse of equation (ay 
comprises only part of the acoustic plan form, the remainder 
heing bounded by so much of the circle 


tery dye ayy tths 


as lies on the wing at time zero. Figure 5 shows the three 
possible acoustic plan forms for points in region 777. The 
limits for the three types are 

Gi) t>r>0 

(ii O>r>- 4M 


(ii) -t M>r>—-1 


camaaeet qe + 
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Piet ee 3. 


Acousde plan forme for region f/f foot thenre o 


and these correspond to the subregions 7//,. /41.. and 111 
identified in figure 4 
potential in, say. region 277, as 


Using equation So. we ean write the 


. 


1 get moe wy tat 
array din : ‘day ‘ 
ory: it ae ap Re 


Sa, 


1 “yeyto "Nog ee 
7 | dy | " “dy i 
oy a ee i Per ter eee 
where 
: M ; : 
Nig oy Line agate OA 


GARDNER'S METHOD OF DESCENT 


Equation (le. governs a four-dimensional s.y.c.f space, 
Our object, of course. is to find for this equation a solution 
that satisfies the boundary conditions in the 2 0) plane as 
specified in equations Zar and ¢2b. Obviously. we can 
always construct a space of more dimensions governed inan 
arbitrary way except that it must satisfy equation (lenin 
an nyc hivperplane. Phen. if a solution tn this higher 
dimensional space which satisfies equations (ar and ibs 
in the say.cé plane ean be found. it represents for & othe addi- 
tional dimension) equal to some constant the solution to our 
problem. This characterizes the method of descent, Tt as 
not obvious. of course, that such oa method lends to any 
simplification; but, with a proper choice of the governing 
equation for the new space, such a possibility. always exists. 

There are examples where various applications of this 
method have proved to be useful. Hadamard’s use of the 
method. mentioned tn the introduction, is chissieat. A simple 
application of his method is the derivation of the velocity 
potential for a source tna two-dimensional stpersante flow 
field. This potential field (which amounts toa step finetion, 
the step occurring at the Maeh waved is easy to derive if one 
field with a dine of sources 
normal to the free stream and uniform in strength. The 
two-dimensional field mentioned above follows minediately 


considers a three-dimensional 


by descent. 

In other examples the additional dimension is measured 
with imaginary numbers and the additional law for the ex- 
tended space is the requirement that the funetional depend- 
ence on the resulting complex variable shall be analytic. 
The method of descending in the latter ease is associated 
with the study of analytie continuation. fn particubar, 
Riesz’s method tdisenussed in ref. 7) for solving equation le) 
Hustrates these concepts. 
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Gardner's method for solving equation Che) is to define a 
five-dinensional space ino whieh a potential function y is 
voverted by the equations 


Wer 
Ve 


ahd show that solutions to equations (0 in this space are 


Wer 
Vu 


ve 
ve 0 


that) 


itha 


general enough to contain general solutions to equation (les 
ina plane & We shall. therefore, proceed by ana- 
Ivaing these equations and eventually let approach a plane 
in whieh the boundary conditions of equations (at and e2h° 
are satistied. For convenience, the latter plane is taker to 
0 plane. 


corstat. 


be the & 
Since equations cas oad (bbe are linear a number of 
possibilities extse for the choice of the dependent variable 
dir oy ot. OL bi Aside from. the 
worn be pur nb where gis the velocity: potential 
of equation Clery: for example, one could det yur ys oO. & 
ee. Gor aca peony OO el ye te be “Phese 
varions choices amount only to relatively minor differences 
If, 


(2a. one 


more obvious choiee 


in the detailed techniqne of the subsequent analysis. 
in impesing the bonndary conditions of equations 
is to use only soureesty pe selutions for both equations (hai 
Therefore, set 


ie Yan. 8] phn re) 
s ic at 


Now differentiate equation (Gai with respect to 
1), 


and itbo. the fast choiee is sufficient. 
(V) 


zand set 


Detining 


Wie ry.t US) 
equation fag ean be expressed in the form 
Wi, W.-W oo i) 
and the boundary conditions for equation (sare given 
directly by equations (20. Phus on the wing 
OW oy so Mt gyy 
oo oer eM (yey 1, { ‘ ) (*) (Liha) 
and off the wing 
ow 
“of gg ALUN 0 (0b) 


Assuming equation (45) to have been solved for the bound- 
ary conditions given by. equations (10), we return to the 
second of the set of partial differential equations (61, spe- 
cifieally. 


Wes $y —¥:7 0 


From equation ¢S), itis seen that the solution to equation 
(9 vields the result 


“Ttean be shown that the solution satisfies the equation 
lini § lim , Uolim _ lim § fim - ' 
so pos [-. racer | e aera tro ee ee ‘ 
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bgt 
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Farther, the boundars conditions for the original problem 
hry oS. fie space require that g be anodd funetion with 
respect tec. and: continuous aeross the O phine except 
Thus go must be zero for. 


evxeept over the wing plan form. 


over the wing phan form, 
The continuation of thas 
condition inte Gry. cg. ft space then implies. according to 


equation 75. that off the wing 


of 
of - 


Henee. both the second partial differential equation and its 
boundary conditions are identier! in form to the first) se 
Hythe, Applying 
to thei dual solution, we obtame the desired 


e veut | Rate 
S fn 


for the potential on a rectangular wing owith oboe b orn 


vive bx equations oY ane respectively, 
equation 7 


result 


supersonie nnsteady motion. 
THE GENERAL EXPRESSION FOR THE POTENTIAL 


The method outlined in the preceeding section wall now 
be applied to obtain integral expressions for the potential 
in any region of the rectangular wine shown in figare 4. 
Consider first equation om) for Wires. tf . Phis equatian is 
the same partial differential equation as that which vowerns 
Further, 
Pespace are identical to those representing a thes 
Since the Miateh 


supersonie steady flow, 
the for. 
planar wing ina steady supersonic flow. 


the boundary values oo 


nimiber ta the steadv-flow analog is 2. 
plan form oof this wing ‘shown tn fig. 6 soa sweptforward 


wing Up baring all supersonic edges Qo e., 


the equivalent 


the component of 

the free-stream: velocity normal to all edges is supersonic: 
Since all edees of the equivalent: wing plain form are 

supersonic. the solution for HW ean be written immediurtels 


Frevere 6.- kquivalent plan form im grt space. 


GRNERALIZED INDICIAL 


th terms of “sources” only. their strength bemg given fy 
fYOae. Thus, ba analogy with the well-hnewn 
results of supersonic wing theory, we have 


bela tious 


; i if wer Mt yada cll, 
ati -— ; 
WIP yb bas 


where ras the areaon the wine ent out by the forecone from 
The analy tie form of WP will dither eon- 
stderably tn oeael of the theee regions above the equisaternt 


Wie. ‘Hy 


Y ea : 


the pow goa. te 


wig shown ain figure 7. 


Prod To Meetous tio whieh anabyoie fort at Wooo dither. 


beeonmes i 
boundary condition for the solation of equation jbo. Thas, 


O plane for whieh yea, g 20, 


The salte or Wo civen ba equation ohh now 


over the portion af the 


: . OF : : 
the varnition of m6 isonow known and tor yO, & 0 
the cosulition 0 apphes UPhese conditions are 
Byres > , 
Al onot sulfierent to determine a umique solution unless 


londing falls to 
> Ngain 
weoobserve that these boundary conditions and the partial 
differential equation cbr are tdentieal to those studied tn 
connection with aostationary planar wing in a supersonic 


the fuether restriction is imposed that the 


zevoons the edge yy Qis approached. toe as ye 


stream As shown tn figure 7. solutions from the t2. € space 
above the &- 0 plane are referred to as Wy. Wy. and W., 
depending on the relation between and 2 in af constant 
Fignre S shows the five different: boundary-value 


Wr. 


plane. 
problems formed by the various combinations of Wy. 


ADA) ? 


FORURS ON 


naan tthe Ae ct Rate ahaa TRL on ee tne a 


WINGS 7 


DEROKMING 
and VE, occurring along constant a hes ii the a. € phine 
and the corcesponding regions in figure 4 for whieh each 
applies. Rach of these five problems as direethy analogous 
to the boundary value problem eneauntercd tn steady-state 
Hfting-surface theory. of a planar. rectangular lifting surtace 
Heo steady supersonie stream. The leading edges of 
these analogous reetangalar plan forms he aloo the foes 


Sf & yr oer se, 4. depending on the value of o 


tovebegeee 


sie tesecareereeeetebe 
Shee ee eee tt sete te. ; 
+o i 
+ 
ssss 
nes pene 
+ 


a er cr cored 


Freere S. The five different boundary-vatie problems mo dia % 


pace, 


ena serene nent eta CL CREE AAIT ate iiss rer theses 
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ad the cade edge hes abeng the dine 7 0. Phetee. tbs 


tueans of this steady flow aiialag, we can immediately write 


the seluten te equatian th. aa the formu 


a Sr TO 


Pe fey 


Wheto od the area of integration a must be discuss. 
Taso possibilities exist for the stipe ofa Furst. uf the 
pert Soy hes te the cht of the dashed dimes ana figure s, 
Whieh ite the aialogois steads-flow problem represent the 
traces of the Miteh cones from the leading-edge tips, ois the 
tromgubu area shows for region 277.0 in tignre 4) part 1 
Hohoweser. soy les between this tine and the sige edue, 
yo Ula asthe trapesoidai ares shown for region Vyo in figure 
Sopart te The hetrer is a well-haown result used in steads 
supersonic ifting-surface theory and first developed bs 
Evin orefos The division of the five kinds of problems 
Ubistrated in figure Sonte the final twelve. represented by 
the reviens i figure 4. es breneht about by the varies 
combinations of Wy Wao and Wy that eam occur in the arena 
asthe port & yoassames all necessary values on the wane 
When y has been determined. the potential in the pli steal 
Plane is fond be equattion To or combining equations «11 


and 2 


| ) fey edy, 
atime © | | ; 
r Joy 


a, a, 


fer 


7 a Mb yode dt 
WJ , 


a 


i (ea? peg te re 
Fite 0 Nrea of ontewration a used in equation p2 


ian ornare a craseamammennasaamnanenammne tomar nana aaaciaa etaciaiiassaaamemantanametimnmmeammematenematatnnata 
gaate ee cement en et nna TITRE Bn weet ee EO fet eS ete ee 


A detuted analysis of equation bs fora pemt a oy. tit 
region Vioof fivire dos wivendn Appendia Ao anda studs of 
this atadysis enables one do write the resalts for abl cet. 
witheut atfhieuttys 


INTER PREPATION OF THE GEESE LTS 


The tesults of the rather aivobved arady sis siven an Qppeenie 
dis Acean be coterpreted ii tertias of the Knew sobutians tos 
sampler boundary conditions  Phese latter solutions tane 
adready been reviewed in oai previous seetion ae whieh of Wats 
shew that the potential ona hifi sarfaee wart ath super - 


wonte edves can he wartter ap the form 


| "por ada iy 
ean || 


From Appendix Actos found Chat the potential ata: geet 
ona reetangubar diftine surfiaee can abwaryss be expressed as 
the sum-of two parts 


ns Ah c rou bit s ay (et ry. 
where 
ee te Ce 
eM, JJ ea ee Ioan 
ae 4 r 
ane 
tore. ; 
Aaa ae wore, Cony das dy 1b 
: ri, d 


The value of roy is given by equation AED in Append 
Acaned the areas of intewrtion. So and Soo are iiastirted: for 
the Saieus regiots Pothronagh Villon fiamre te 

Letts first inspect equations 15 an light of their posstble 
analogs with the familiar solution for the steady-state. ree. 
tanvular lifting surfaee  Efoa reetanpular wai diaaving ar- 
bitrars twist and: camber as placed inca steady sapersente 
flow. the solutian foe the potential on its surfer earn afse fe 
expressed as the sum of two parts 


cir ay ne a on ee Tey? 
where, if 
rr a coe 
1 | 5 wi da dy, - 
cna : yout 
ri, , 
and 
PCC wheat z 
ey || ‘ " cbTb 
: Li 7 


These equations ean be constried in the follow ser siapde wai 
Equation [Tas represents the potential indeed ats 0 ba 
datribation of sonrees over the wing plate forte. eaek soarce 
having a streneth proportional te the loca) streamwise slope 
of the Upper surface Phe area Soo as shewr on Here 10 
the portion of the wang wither the Miteh forecane fron 
sy) Equation tb dias a simile iaterpretation ot aise 
represents a distribution of sontrees aver the ware ened 
having a strength proportionai te the locab slope ot tlie upper 
surface Bat the area of integration Sos sew that portion 
of the wee within the Misch forecone from the part. ga: 


GENEKALIZED INDICEAL FORCES ON DEFORMING WINGS y 


bietae Lt Neoustie phon form: for port in recite Ve fanare 


“Hinton identical to that existing between! and Sou the 


steadvestate case. oti other words, S, ts the acoustic plan 
form for the event ao,0.0 and Sots the acoustie plan form for 
the event a. yt Phasteally So represents the portion 


of the wing's lower surface containnie disturbances whieh eat 
Prove TO Aree Sted S cised un eaqumtnors TT 


; , : surface At this point the similarity between the stead ane 
Chats. within the cone which forms a mirror image of the : 


ane unsteady solutions ends since Che iifhtonee of the slopes ta 
plivsterl Mitel foreconme in the vertteal phine containing the 
wings stde cedure The potential e¢ 2a0) represents the 
HHifferenice between the potentials fora wing with a verticoally 
sVinmietrieab thiekness distribution and oa surface with ne 
Hiekness having the same starpe as the upper sarfiee of the 
nordiftinee Wittig, 

Lets ceturn now te equations Toe. dusts in the steady - 
state conse, oa OQ represents the potential induced at 
nay by occbistribation of sourees cee eq. over the wing 
phir fort. enech proportional tar the level slope of the WIT, 
Hat tow. since the wing ts in metion, with the added eou- 
dition that they be locrl slopes at the appropriite time. The 
aren S.. shown a figure PEO as just the acoustte phan fern 
defined earlier in the discussion of equations Si and «4 
Physteally S, represents those pomts on the wing from whieh 
chsturbanees ean. at the time A inflirence the flaw at suo. ft 
tothe generalization. in the stationary coordinate avstem, of 
the wing area boanded by the Maeh forecone. 

Pipa Ot: and 2? cry0.f is similar 
to that between there steady-state analogs. “Phus. again, 


The relation between ¢ 


2 orth represents the difference between the potentials 
foram uneambered nonlifting wing and oa difting surface 
having the same shape as the top of the nontifting wing. A 
more striking similarity lies in the relation between S, and S,. 

We have already seen that S, is the acoustic plan form, 
andoasat turns out. Ss the reflection of the acoustic plan form 


in the cortical plans containing the side sdqe tsoe fig. 12) a figure 4 


wat the time © tutluence the How at voz00on the wines upper 


Fyectre 12 Reflected seaustio plan fort for pita im regia Way 
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the reflected phan fort is not the same as it is for the slopes 
in the basie acoustic plan form: the influence in the former 
case bow being given by the integral Curj.y.) detined in 
equation CATO), 

One can show, by simply referring the results given in 
equations (13) tea coordinate system fixed on the wing, that 
equations (Man oand Clb) are identical, respectively. to 
eqtationus cl7ab and ob) when they apply to regions UZ 
and V7 in figure ¢: regions in whieh, for indicial-type motions, 
the flow is steady relative to the wing. 
Chontand cobs extend Evvard’s vreflected area” concept to 
all partsof a rectangular wing tn supersonic unsteady motion! 


Henee, equations 


THE GENERALIZED FORCES 
REVIEW OF LAGRANGE’S EQUATIONS OF MOTION 


In order to define more clearly the subsequent: concepts 
and notation, we will brietly review Lagrange’s equations of 
motion as applied to distorting wings and will examine a 
simple application toa reetanular wing. 

Lavrange’s equations are usually written 


dot of ol’ 
I oq. aq." Oy. a a a cs 
where 

TO kinetic energy of the wing 

f° potential energy of wing 

2 aveneralized cexterual: force 

yd generatized coordinate 

Tnthe present appheation q is the amplitude ata given time 
ofa polynomial measuring 2. the vertical displacement of the 
wing's camber tine from the > Oplane. Phus. relative toan 
voi, Coordinate system that is fixed on the wing, see figure 13 


her yt Sy to agg (1a) 


s 


+3 


s 


= 
fas 


a3 


Breve 13. 0 Wing in moving coordinate sv-tem, 


Sdrosof farther interest to notice that equation Chats cin be reduced toa double mteser 
invelw ing warts by using, for example, the transformations ¢-1;© Wrrand r=/—f amd 
‘otegrating With respect tag, 
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The wing's Kinetic energy can be written 


age ea, 
7 ff honey daddy (20) 
< 
Where a ots the wing nitss per unit plan foria area. Using 


egnation (fs we find 


: pe 4. | | Por ya cyyiner yd ady 

AO De oe 
ol ei 
OW- 


The potential energy is usually dithicul te evaluate analyst. 
ieally. However, it can often be determined experimentally 
das Will be seen: ba measuring the frequencies af the free 
For the present assume Chat the wang isa 
The potentnel 


Vibration modes. 


homogencous plate of constant thickness. 


energy for sucha wing can be expressed as cref o 


. oe eres Oh oh oh 
{ > {| {« Ayo ot oy & “as ~( Sron ) }} da dy, 


(24 
which leads to the equation 
ole "Te. . Lor. ol’. 
NY sp > » ; 
a Dy 4 If [save on (., Sos 


lola? oF ol’. ) Water & 
2Or- Oye O04, OF OY va 


te 


where yois Poisson's ratio, Te Of Ore - Of Ove. and 


p 7 Younes modulus) (plate Ghieknesss 
, fet 
Sib we) 

Now. if the ceneralized coordinates have been portmalized 
so thantench measures the amplitide of a free vibration mode, 
all terms in equations ‘2h) and (23 
the product of 2. and 2 oare zere. 
such normalization. we ean write 


i-| | Pour yamayodaedy.- Dy, | | {irr 


OP ONe pO TY gig to ns 
as bee on aroy,) | peat > te Vee 


(24s 


Involving the integral of 
Assuming, henceforth, 


Finally. dividing through by the coefficient ef 7. and ex- 
pressing a generalized foree as the integral over the wing plan 
form of the product of the rth mode shape and the loadings * 
Spe. induced on the wing ba each of the mode shipes 
considered, we find 


qu. SS Pasge ( Ba ) dary. 


s ros 


SS Por yomian ap idarydy, 


de Qeer 


where a, is the frequeney of the cth free vibration mode. 

We will write sap . Phas os 
possible withown a confiisien of notattet since the ge 
apd excdiude Phe term ye. 


Where qos the fr 


YAP pn cstreau da naanie pressare 
rl zed eoortinates ape expressed as 


BA 


GENERALIZED ENDICIAL 


If the free-mode frequencies are experimentally deter- 
mined, equations giving the wing's 
potential energy, never have to be evaluated, Further, in 
such cases, equation 123) applies to quite general wing struc- 
tures with varying density. Usually in the applieation of 
equation 25), one uses the actual frequency w, of the free 
mode but, in evaluating the aerodynamic 
analytical expression that only approximates the rth mode 
shape. Let us examine the eneralized force tert in equa- 


tion (23). taking. for simplicity. only one term of the sum: 
£ , A 


such as equation (23) 


forces, uses an 


ay ) dryly, (20) 


qo 


QQ: uw» gee ( 


According to what has gone before, the mode shape poly- 
nomial Piya) has the form 


while (Ap qu), is the loading cocthcient corresponding to an 
indicia! deflection (see previous section on boundary condi- 


tiats) 
oe Pad WETS GAY Ee Mt 3 
h 1. HL) ( =) [( ¢ ) +1( 2% c ) (28) 
which gives a vertical velocity distribution 
: Ys » 
we’ gel) (= ) (BY (29) 


Now a generalized indicial force coefficient can be defined as 
follows: 


Sere 1 ny (usy 4 
Wait) = cael » ff( AY (2) tap qededdevtn 80) 


(The calculation of these quantities £17(0 will be elaborated 
in the Sines the generalized force (, 
intended to apply to any motion, not necessarily indicial, 
it is necessary to apply Duhamel’s integral to the indicial 
force coefficient M7. thus, 


d ak od rer! ) 
GO, ys qd zr’ dy’ 
BRS Ede Je ql) 

As ain example, consider now a simple one degree of free- 
dom vibrating plate. The plate is fixed to a wall and 
restrained along its leading edge. The mode shape is as- 
sumed to have the form 


h- eq) (2 2) ( Sy, (32) 


so fora plate with uniferm density and thickness 


*o we mse ps\i 
a2 = ed 
- J -8 asf teat ‘; (454s) 25 ( C ) 


Equation (25) now becomes 


next section.) 


(31) 


ey (33) 


ian 8 
qr orn 


FORCES ON 


WINGS 11 


DE FORMING 


For this case. we have the generalized indtcul foree coefficients 


TU, and FEU, 
aoe 
( ye Pail! 
Hey J od ee 
yz: ftir’) 
Ryn ae oo {7 
ee re iee 4) 
t 4 Ip 
Therefore, equation (337 ean be written 
aa -. maa TT i 
yi wr i i 
“Aor 
et r 
U cT (33) 


THE GENERALIZED INDICIAL FORCE COEFFICIENT 


Lt is clear from the previous section that a study of the 
dynamic behavior of rectangular wings moving at supersonic 
speeds can be carried out if ene can obtain values of the 
generalized force coetticient, fact: as defined by equation 
(0). We will now show how these values can be obtained 
from the solution the aerodynamic boundary-value 
problem represented by equation C14}. 

It was convenient in developing equation (14) to use a 
coordinate svstent- y.c.6 whieh was fixed in space so that 
the left edge of the wing moved along the s axis as shown in 
figure t. On the other hand, in’ studying the dynamic 
problem it was more convenient to Use an 43.45.2306 system 
which is fixed on the wing. In order to convert the results 
in one coordinate set to the other, let us first transfer results 
in the say.c.Oset to the 24.44.74. 0set (shown in figure 14) and 
then, finally. Cransfer to s3,%.2u¢ coordinates. 


to 


Lo t4tM 
Vs 74 


23724 
tet 


44 


Fietre 14.—-Transformations from moving to fixed coordinate system, 


i t ne 
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Phe mdicial force coettcient Ais is detined as follows: 


Par am NO taal Be te ar 
Fey ad val dy (7 yc oe sth 


In order to transfer the axes from the set shown in figure 1 
to the more conventent set of figure 14, so that mode shapes 
are svminetric or usvmmetric about) the wing's spanwise 
center fine and the force coefficients denoted (44 can be deter- 
mined, we proceed ax follows. First, the loading coeificient 
fora wing in the oye system with dewnwash given by 


wy por Mty'yy-s 
es -( : ‘ye ) 
eee GY CY 


uel < 


” 


is obtained, 
sum: 


This loading coefficient: can be written as a 
in 


P Nu ai \"" he 
(Go) me Bol yay) 


wad 


Now the quantity {37 is defined in the xy. y.,2a¢ system as 


Bone Sg eS te OEM AY 
ft ae s ( i: Cx )( ys yi = )" 


This fast integral can be written as 


Nae fy OMY * ay 


Qo 


eB HVAT HEC) 
“ eae de fa (7) Mty! (Yy(22)" 


By using equation (36) we find 


fips ae "| Ss 


Sra) zB 
()(a) Fe on 


where all forces are responses to a unit indicial disturbance. 
Note that if equation (37) is applied in the case of a wing 
cantilevered on a wall, both » and g must be even in order 
to satisfy the boundary conditions of reflection in the wall. 

By superimposing boundary conditions and their result- 
ing solutions, one can further show that the value of 2 
given by equation (37) is valid for all reduced aspect ratios 
3.1 greater than 1 in spite of the fact that the value of Fir 
given by equation (36), as it stands, applies only to wings 
for which 3-4 is vreater than 2 

Given H7), one can determine the generalized force as- 
sociated with the generalized coordinate q, by means of the 
superposition integral as illustrated by equation (34). 


ALS ISORY 
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DETAILS OF CALCULATION 


The details of actually evaluating the indicial force co- 
effictonts from the solution for the potential presented in the 
first part of this report are discussed in Appendix B. 
siderable 


Con- 
labor is invelved in) such calculations, and an 
attempt was made to discover recursion formulas by meas 
of whieh certain derivatives, for the reetaugzuher wing, could 
be expressed as combinations of others. ‘This attempt was 
successful and vielded the following results: 

Consider equation (360. Tategrate the a integral in this 
equation by parts, setting 


*. in 
url | jer dy: deve) tae Miyata 


Yu 


Then, since by equation (B7) in Appendix B 


O ap fap il" cy 
Or qu c Yo ’ ‘ 
one finds 
5 ! y ‘ n + * 
Fy aes | ee tthe ae (38a) 


Inspection of equation (37) shows that the same relation 
holds for the generalized indicial force coefficients (44; that is, 


(38h) 


From this relation, it is seen that only the forces #4" need 
be determined by integration; the forces for higher values of 
the index J can be found by combination of results for dif- 
ferent values of the mode shape index /. 

As a simple illustration of the results presented so far, we 
can calculate the indicial force derivative for the cases |= 
n=g=0,j)=0, 1. The case ;=0 corresponds to the indicial 
lift coefficient for a flat, sinking, rectangular wing, and the 
case for j=1 corresponds to the indicial pitching-moment 
coefficient for the same wing. Sinee n= g=0, equation (37) 
gives 


oe = FM, 


Thus. with j=0 and identifying ao Oo as angle of attack 
a, one finds from Appendix B 
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Next with job. and using ©,’ to designate the pitching 


thoment measured about the leading edge of the wing, 


23 i ae on 
‘ me =( ape t i ) fs TM {( I~, tn }s 
yy ne pnt} OSS yp ' 


2 fi t; oy Mta~ 1 
Tae [( ae J eos t " + 


t + Mt, yt ~ Mtr’ | 
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ora ee 
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38 
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ie ‘ | 
=a saa ME yp 
These expressions agree with those given by Miles in refer- 
enee m, 
The above results can be used to demonstrate the useful- 
ness of equation (38a). ‘Paking g=n=q=0, (==) in that 
equation gives 


Fin= Fyo— FX 


or. for the present case, 


which represents the equality 


14 = +O, 
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Combining, we find 
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20 =t 5 (— Mt) 32 A sae+syt OStS 74, 


P set 1 3 : pot 


13 


WINGS 


that is, the lift coefficient fora pitching wing equals the sum 
of the lift and pitching-moment coefficients of a sinking wing 
(primes indicate the wing is pitching about and moments are 


Heuer, 


measured about the wing leading edges. 


i MMt. 1 MM ; 
af “ (1+, } cos I, + 3 cos "CM - Bhat 


SMM oe ayant Ui eet ea ei 
2 yir-- a gene BEML Dts + 


2 ty wl | | 
ar -rin] be we ShS apy 


2h yo 4 a piadl 
“p sp. fo “> M-A-1 


A further application of equation (38a) provides the piteh- 
ing-moment coefficient for a pitching flat rectangular wing. 
Thas, with /= j= 1. n=4=0, equation (38a) gives 

n= CRA Pa 
which becomes 
A= FO — fe) 


and so 
=e j a) 


From equation (B21) in Appendix B it is found that 
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or relat r » generalized indici toes fini U = i ‘ ; LX . il l 
Another re lation among the generalized indicial forces.’ Si bee 1 ( Ne, ‘ep 1° Ay2G 1>( ) fis 
w j=0 ry 


ean be derived by means of the reciprocity relations given in rea J 
reference 5.) The details of the derivation are given in 
Appendix C and there results 


wet t= nel 


140) 


The inner sum on the left ean be evaluated. Thus one has 


i 
I 
| 
1 
l | 
Jv . 1 
Sew yan (— Le ( ys (39) 5 v 
Pex) y uso en | we [be tha? Sie (P)a- ve 
“ 
Equation (39) can be used in two ways; one, as a means | 
for checking the internal consistency of a set of calculated | 
. + eos { 
generalized indicial forces, and the other, as a means for | 
expressing a given force in terms of a set of others. 
Conskler, as an example of the former use, the case for 
which Jej7- 0. Then 


| rai wer 


Equating cocthicients of 7, 


ass ' 
Pay An | 
. : A a | fo 0; rep 
From equation (37) we can express this relation in terms of | Sot i» ( *)(*) Cpe 
the calealated quantities F'" thus oe . wee 
: 
, caer and equation (40) becomes 
See (Y Baw (YP ee | 
oe 0 : 
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Bie (Z) ecm CYS Fy 
rau a= = CONCLUDING REMARKS 


If now note 3 the following relation results A method is presented for evaluating the generalized forces 


ona rectangular wing flying at supersonic speeds and having 


ys Mdy if pr . . . 
(Fv W+5 eae Fy) + 3CF 33 Fo) + an aspect ratio such that 8421. The generalized coordi- 
nates used to define the wing’s behavior are the amplitudes 
aS 9) (FR — F%) +2 (3 yA mF) of downwash distributions expressed in terms of polynomials 
" inv and y, the chordwise and spanwise directions, respec- 

tively. 


whieh provides a useful check on the computed quantities. 
Next let us solve equation (39) fora given foree. Perform 
the sum operation 


Numerical results are presented in table [for generalized 
indicial forces on a wing having an aspect ratio of 4 and 
flying ata Mach number equal to 1.1 and 1.2; the polynomial 
a 7 coverage being O</< 1 and O<n <5, where ww rly", 
> (=n ( ) 
j=0 

Ames ABRONACTICAL LABORATORY 
on both sides of equation (39), and reverse the order of NaTIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
summation on the left side. There results Morrerr Fretp, Cacrp., June 30. 1954 
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APPENDIX A 


EXPRESSIONS FOR THE POTENTIAL 


In order to write the expressions for the potential in all 
regions shown in figure 4, it is sufficient to derive in detail 
only that for region Vo Having carried out this analysis, 
one ean determine the expressions for potential in’ other 
regions without difficulty, 

Consider, therefore, equation (13) and let o and 7 apply 
to region V,.— First, it ts necessary to determine the poten- 
tials Wy and WW’, in the (a2,é space. From equation (11), in 


conjunction with figure 7, it is found that 
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| With the values of Wo given in equations (ALi and GA2) 
| itis possible now to solve equation (6b) for y. figure 8 giving 
| the required data in the &y plane. Thus, if Rf (gg) 
hy -me 
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Now apply the operation of equation (7) and the potential gy, is given by 
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where Ry? =§?—(y—y,)? and the bars on the integrals signify 
that the finite part of the integral is to be taken in the sense 
defined! in reference 10 and that the order of integration 
cannot, in general, be reversed.? 


For convenience set 


al, (A5) 


Ww 
Y= 
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1 
5 
where J, is the nth integral group on the right-hand side of 
equation (A4), 
Consider the first of these integral sets. Using equation 

a 

1 For the subsequent analysis to hold, the definition of the finite part given in reference 10 
is essential. This definition differs from that given by Hadamard when it applies to muitiple 
integrals, 

? Since the order of integration plays an important role in the following development, 
integration first with respect to r and then with respect to y will be denoted Sdy S dz f(r, ») 


while integration first with respect to y and then with respect to r will be denoted fdr Sdy 
ty). When the notation SS f(z, y) dydr is used, the order of Integration is immaterial. 
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In order to simplify this expression, the order of these in- 
tegrals will be rearranged so the integration with respect to 
£, can be carried out first. Phe technique of changing the 
order of repeated integrals with strong singularities set forth 
in reference 10 will be used here. Consider the change of 
order in the &,2; plane. Pretend for the moment, that the 
t, integration has been carried out. Then the highest order 
singularity (since w, is bounded) in the &, +; plane has the 
order 3/2 which is weak in the sense that no residual occurs 
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when the sequence of integration is reversed. The top of 


tigure 15 shows the area of integration, so immediately 
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Fietre 15.—Areas of integration used in analysis. | 
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figure 15. Tn this case an inherent singularity exists at the 
confluence of the singularity lines of the qitegrand > mamely, 


where &, You und ¢, + yin at ES. 
order can therefore not be performed directly, bat account 
must be taken of the existence of a residial term csce ref. 10. 
This vesidual is defined as the difference between the two 
integrals taken in different orders over a vanishingly simall 
tthe dotted 


The residual 722, is then, 


The change of 


region surrounding the inherent singularity 


region in bottom of figure 13). 
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where ré> (r—20,)> (y-yyt. The second integral vanishes 
(see ref. 10), and, passing to the limite. +O in the first integral 
there results 
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where the square brackets again mean that the retarded 
value is to be taken. Thus, the integral J, can be reduced to 
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In the same way, the integral /, can be reduced, and 
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which is recognized as Kirchhoff’s formula, equation (3), 
with an acoustic plan form bounded by the cirele 
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The reduction of the integrals /,. 4. /;, and J, is quite 
similar, leading to the sum 
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Examination of the limits on these integrals shows their 

total area of integration is that shown in figure 11. But 

this area corresponds exactly to the acoustic plan form S, 

for a point in region Vi! Hence, denoting the combination 
of terms in equation (A7) by ¢" we can write simply 
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It now remains to caleulate the integrals /; through Jj. 
Designating their total effeet on the potential, ¢? one can 
readily show (since no inherent singularities arise in’ these 
cages) that 
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GENERALIZED INDICIAL FORCES UN DEFORMING 
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In terms of this expression, equation (AQ) ean be written 
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where the area (S.)y, is illustrated in figure 12. 

In order to give expressions for the potential in every 
region of the wing shown in figure 4, one can show that it is 
only necessary to vary the areas over which the double 
integration in equations (A8) and (AIL) are carried out. 
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and only the acousnc plan form 3, changes with the region. 
In the case of ¢ 7. the part of the potential due to the exist- 
ence of the side edge of the wing. equation (ALD) can be 
generalized and written 
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where the integrands are defined in every case by equation 
(AlQ) and only the “reflected” acoustie plan forme S, 
changes with the region. The region S, is always bounded 
by portions of the “reflected” cirele 


irre dyer ee E 


and the “reflected” ellipse 
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This is evident in connection with the source portion yg”. 
for in every case A ; ; 
1 [wl Figure 16 shows sketches of both S, and S, for all regions in 
e ~~ ff ro dryly, (A12) figure 4. The absence of a sketch indicates that the corre- 
Se sponding integral does not exist for that region. 
APPENDIX B 


THE GENERALIZED INDICIAL FORCES 


THE LOADING COEFFICIENT 


In order to determine total forces acting on the wing, it is 
tirst necessary to obtain expressions for the loading coefficient 


Ap qo. According to the linear theory 
Ap_ 4 ow 
qo CoM Of en 


so it is necessary to differentiate each of the expressions for 
potential. As an example, consider, as in Appendix A, just 
region V, of figure 4. The loading coefficient will be divided 
into two parts Ap? /q and Ap /qo to correspond to the po- 
tentials ¢ and ¢®. Thus, using equation (Al1) 
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since the derivative passes the 2).y, integration without 
effect. Referring to equation (A10) for the function Cr, 1) 
we next find its derivative with respect to ¢. 
then for 4,<0 
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Notice that if w, does not depend on (7,4 A/t)) the integral 
term in equation (B3) vanishes, while if it does, then the 
integrated term is zero. Next, for 2, >0. 
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In this ense, both terms exist: unless aw, ts not a funetion of 
(2, + Mt,), in which case the integral vanishes. 

Substitution of equations (B3) and (B4) into equation 
(B2) will now vield an expression for the loading coefficient 
corresponding to the influence of the side edge; 
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Firetree 16.--Sketches of areas of integration, S. and S,, for all regions Ficgvre 16° Continued 
in figure 4. 
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GENERALIZED INDICIAL FORCES ON DEFORMING WINGS 1Y 
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Its clowr that. even for small values of the indices J and a, 
the required integrations for the determination of total forces 
on the wing pose formidable problems. There ts, however, 
a property of the loading coetheient corresponding to vertical 
velocity distributions of the type chosen here ceq. 2h) that 


will iatectally shorten the requisite labor. This may be 


expressed as follows, adopting the convention that Ap!” qo 
corresponds toa downwash distribution proportional to 
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DETAILS OF EVALUATING THE GENERALIZED INDICIAL FORCES 
In caleulating the generalized indicial forees by means of 
equation (36), i has been shown that only the value zero 
need be taken for the uidex 4. Thus we must find 
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The values of the loading coetticient Spe" ga are found by 

differentiating the expressions for potential given in the test 

bart of this appendix, 


[tis convenient. in evaluating equation «Bao to eousider 
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Hos found that 2 seems to bave different: representations 
according to the interval in whieh dies. ‘These expressions 
can, however. alb be expressed by the same formula. The 
portions of ZL corresponding to the parts eo and ge? oof the 
potential are similarly siguitied. and we have 
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wand [a 2] means the greatest integer contained ina 2. The 
function Jig) may be expressed as summations, and it his 
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Where ( c is the binomial coefficient 
CP iis: 
277 Mp2)! 
and Bop. qi is the beta funetion 
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The funetion Ji. 0) has been calculated for gq. on taken 


O, 1.2.8. 4.5. Beeause of the property (B13), it is only 
necessary to give a triangular array, whieh appears in the 
above table. 

Now consider the functions Aya) and Ayiv), defined after 
equation (BI2). Tt is convenient, for computational pur- 
poses, to express these in terms of the incomplete beta 
funetions, defined as 
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A tabulation of the incomplete beta funetions is available in 
reference V1. Note that when the symbol 2 is written 
without a subscript, the complete integral is meant. that is, 
in eqtation (BiSbi. 2 equals QO. Tt is necessary to exercise 
some care when interpreting Auiy) and Aytvt as) beta 
funetions because of the upper limit. Thus, sinee 
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The integrals Jivs and /yivrean be simplified by reversing: 
the order of integration. This can he accomplished ino a and, sinularly, i ean be showre that 
strught-forward manner by merely inspeeting the region of: ea 
Integration in the su4 plane. Consider first the integral: Pi ix M 1 | Mat Soeraiie 
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DERIVATION OF RECIPROCITY RELATIONS 


According to reference 5, the reciprocity relation for 
general three-dimensional unsteady motion can be written 
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